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1. INTRODUCTION 



'g\ ' Gluonia, Scalar and Hybrid Mesons in QCD* 

^ ! Stephan Narison'^ 

Q_|! ^ Laboratoire de Physique Mathematique, Universite de Montpellier II Place Eugene 

^ I Bataillon, 34095 - Montpellier Cedex 05, France. E-mail: narison@lpm.univ-montp2.fr 

CN I For some experimental guidelines of the next millenium, I review the determinations 
^ ' of the masses, decays and mixings of the gluonia, scalar and hybrid mesons from QCD 
,—1 ■ spectral sum rules and low-energy theorems, and compare them with the lattice. 

> : 
o . 
i> ■ 

^ : 

0^ , Since the discovery of QCD, it has been emphasized that exotic mesons beyond 

^ I the standard octet, exist as a consequence of the non-perturbative aspects of quantum 

0> ' chromodynamics (QCD). Since the understanding of the nature of the r]' [^, a large 

r-| ■ amount of theoretical efforts have been furnished in the past and pursued at present for 

5^. predicting the spectra of the exotics using different QCD-like models such as the flux tube 

Qh! , the bags [Q, the quark [|| and constituent gluon models . In this talk, I shall 

^ I review the present status of the predictions from the QCD spectral sum rules (QSSR) a 

■ la SVZ [0] (for a review, see e.g.: |P) and from some low-energy theorems based on Ward 

' identities, which I will compare with the lattice results. 

.c^ ; 2. QCD SPECTRAL SUM RULES (QSSR) 

2.1. Description of the method 

Since its discovery in 79, QSSR has proved to be a powerful method for understanding 
the hadronic properties in terms of the fundamental QCD parameters such as the QCD 
coupling as, the (running) quark masses and the quark and/or gluon QCD vacuum con- 
densates. The description of the method has been often discussed in the literature, where 
a pedagogical introduction can be, for instance, found in the book In practice (like 
also the lattice), one starts the analysis from the two-point correlator: 

i^Hiq'') =ijd^x e*^^' (0|TJh(x) (J^/(0))^ |0) , (1) 

built from the hadronic local currents Jh{x), which select some specific quantum numbers. 
However, unlike the lattice which evaluates the correlator in the Minkowski space-time, 
one exploits, in the sum rule approaches, the analyticity property of the correlator which 
obeys the well-known Kallen-Lehmann dispersion relation: 

^H{q') = r , t • - Im^H(t) + (2) 
JO t — q'^ — le TT 
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where ... represent subtraction points, which are polynomials in the g^-variable. In this 
way, the sum rule expresses in a clear way the duality between the integral involving the 
spectral function Imipuit) (which can be measured experimentally), and the full correlator 
ipnil'^)- The latter can be calculated directly in the QCD Euclidean space-time using 
perturbation theory (provided that — + {m being the quark mass) is much greater 
than A^), and the Wilson expansion in terms of the increasing dimensions of the quark 
and/or gluon condensates which simulate the non-perturbative effects of QCD. 

2.2. Beyond the usual SVZ expansion 

Using the Operator Product Expansion (OPE) 0, the two-point correlator reads: 

Mq')^ E ^ C{q',u){0{u)) , (3) 

D=0,2,4,... [ — Q ) dimO=D 

where u is an arbitrary scale that separates the long- and short-distance dynamics; C are 
the Wilson coefficients calculable in perturbative QCD by means of Feynman diagrams 
techniques; {0{iy)) are the quark and/or gluon condensates of dimension D. In this paper, 
we work in the massless quark limit. Then, one may expect the absence of the terms of 
dimension 2 due to gauge invariance. However, it has been emphasized recently that 
the resummation of the large order terms of the perturbative series, and the effects of the 
higher dimension condensates due e.g. to instantons, can be mimiced by the effect of a 
tachyonic gluon mass, which might be understood from the short distance linear part of the 
QCD potential. The strength of this short distance mass has been estimated from the e~^e~ 
data to be ~ -(0.06 ~ 0.07) GeV^ which leads to the value of the square of 



the (short distance) string tension: a ~ — la^A^ ~ [(400 ± 20) MeV]^ in an (unexpected) 
good agreement with the lattice result |jl2| of about [(440 ± 38) MeV]^. The strengths of 



the vacuum condensates having dimensions D < Q are also under good control: namely 
2m{qq) = -mlfl from pion PCAC, (a.G^^ = (0.07 ± 0.01) GeV^ from e+e" J = 1 
data |lO| and from the heavy quark mass-splittings [jl3|, Ois{qq)'^ — 5.8 x 10~^ GeV^ [p!o|] , 
and g^{G^) ^1.2 GeV"^ {asG^) from dilute gaz instantons 



2.3. Spectral function 

The spectral function is often parametrized using the "naive" duality ansatz: 

- ImiPnit) ~ 2M]pf^6{t - Ml) + "QCD continuum" x ^(t - tc) , (4) 

TT 

which has been tested [^ using e'^e^ , r-decay data, to give a good description of the 
spectral integral in the sum rule analysis; fn (analogue to /^r) is the the hadron's coupling 
to the current ; 2n is the dimension of the correlator, while t^. is the QCD continuum's 
threshold. 

2.4. Form of the sum rules and optimization procedure 

Among the different sum rules discussed in the literature |^ (Finite Energy Sum rule 



(FESR) |T5[, r-like sum rules |T6[,...), we shall mainly be concerned here with: 



• The exponential Laplace unsubtracted sum rule (USR) and its ratio: 

Cnir) = dtr exp(-tr) - Im^nit) , 7^„ = -— log£„ , (n > 0) ; (5) 

JO TT ar 
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• The subtracted sum rule (SSR): 

C-iir) = / - exp(-tr) - Im^jf(t) + ^^(0) . (6) 

Jo t TC 

The advantage of the Laplace sum rules with respect to the previous dispersion relation 
is the presence of the exponential weight factor, which enhances the contribution of the 
lowest resonance and low-energy region accessible experimentally. For the QCD side, 
this procedure has eliminated the ambiguity carried by subtraction constants, arbitrary 
polynomial in g^, and has improved the convergence of the OPE by the presence of the 
factorial dumping factor for each condensates of given dimensions. The ratio of the sum 
rules is a useful quantity to work with, in the determination of the resonance mass, as it 
is equal to the meson mass squared, in the usual duality ansatz parametrization. As one 
can notice, there are "a priori" two free external parameters (r, tc) in the analysis. The 
optimized result will be (in principle) insensitive to their variations. In some cases, the 
tc-stability is not reached due to the too naive parametrization of the spectral function. 
One can either fixed the tc-values by the help of FESR (local duality) or improve the 
parametrization of the spectral function by introducing threshold effects fixed by chiral 
perturbation theory, in order to restore the tc-stability of the results. The results 
discussed below satisfy these stability criteria. 



3. UNMIXED GLUONIA CURRENTS, MASSES AND COUPLINGS 
3.1. The currents 

In this paper, we shall consider the lowest- dimension gluonic currents that can be built 
from the gluon fields C^^ and which are gauge- invariant: 

u,d,s 

a 

'IT 



d,A^{x) = tr + ^ m,g(^75)g , Js = gfatcGl^G'^^^G'^^ , (7) 



u,d,s 



where the sum over colour is understood; 6'^ is the trace of the energy-momentum tensor 
Ofiu] dfj,Af^{x) is the U{1)a anomaly; rUg is the light quark mass and /3 is the /3-function. 

3.2. Masses and couplings 

The unmixed gluonia masses from the unsubtracted QCD Spectral Sum Rules (USR) 
r?! are compared in Table 1 with the ones from the lattice |T^p!H| in the quenched approx- 



imation, where we use the conservative guessed estimate of about 15% for the different 
lattice systematic errors (separation of the lowest ground states from the radial excita- 
tions, which are expected to be nearby as indicated by the sum rule analysis; discreti- 
sation; quenched approximation,...). One can notice an excellent agreement between the 
USR and the lattice quenched results, which the mass hierarchy: Mq++ < Mq-+ ^ M2++ , 
expected from some QCD inequalities However, this is not the whole story ! 



4. PSEUDOSCALAR GLUONIA 
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Table 1 



Unmixed gluonia masses and couplings from QSSR [0 compared with the lattice. 





Name 


Mass [GeV 










VT, [GeV] 


fn [MeV] 






Estimate 


Upper Bound 


Lattice 


12 








0++ 


G 


1.5 ±0.2 


2.16 ±0.22 


1.60 ±0.16 


2.1 


390 ± 145 




3G 


3.1 


3.7 








3.4 


62 


2++ 


T 


2.0 ±0.1 


2.7 ±0.4 


2.26 ±0.22 


2.2 


80± 14 


0-+ 


P 


2.05 ±0.19 


2.34 ±0.42 


2.19 ±0.32 


2.2 


8-17 



4.1. Testing the nature of the E/l 

We test the gluonic nature of the E/l by determining its decay constant /t, from a 
saturation of the USR (Eq.||) and SSR (Eq.D pseusdoscalar sum rules by the rj' , E/l and 
the gluonium P. One obtains l]!^ a value of consistent with zero, in agreement with the 
estimate from the J /if) ^ ''y + E / l decays (/^ ~ 7 MeV), and smaller than k, 30 MeV 
20[] and /p ~ (8 ~ 17) MeV. The quarkonium-gluonium mass mixing angle is determined 



to be small [Op 12°), from the off-diagonal two-point correlator [pT] ,p| (see also ||22|| ): 
^,,(g2) ^^jd^x e'^^0\r tr G„^G"^(x) ^ m,q{t^,)q (0)|0) , (8) 

u,d,s 



which, then, suggests a small mass shift of the physical states after mixing. 
We may conclude that the E/l is likely the radial excitation of the t] or/and v( . 

4.2. Radiative decays of the pseudoscalar gluonium P 

Using \Qp\ ^ 12°, and the OZI rule, one can predict pT|J^,p^: 

r(P ^ 77) ^ (1.3±0.1) keV , r(P ^ p7) ^ (.3 ± 0.1) keV , (9) 

where the errors are probably underestimated. This result is testable at BES. 

5. TENSOR GLUONIA 
5.1. Spectrum 

From Table 1, the lowest ground state and the radial excitations of about i/t^ are almost 
degenerated in masses, which suggests a rich population of the 2"'""'" gluonia around 2 GeV. 



Though the C(2-2) is a good gluonium candidate [^, it may not be the lowest ground 
state. A complete analysis needs systematic scannings of the region above 1.9 GeV and 
further tests of the old BNL candidates qt- 

5.2. Quarkonium-gluonium mass mixing 



A QSSR analysis of the tensor correlator p^j8[], similar to the one in Eq.(H), leads to a 



small quarkonium-gluonium mass mixing angle of about —10°. 

5.3. Tensor T decays 

One starts from the universality of the vertex form factor 



(7r(p)|6'^i,|7r(p')) ~ (Lorentz structure) x 1, at g = (p — p') = 0, (10) 
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and write a dispersion relation. Using the /2 nn data, one can deduce 



r(r ^ TTTT + KK + rjrj) < (119 ± 36) MeV) , T{T -> txtx) ^ 10 MeV < 70 MeV , (11) 
in agreement with present data [^. A non-relativistic relation between the 0+'^ and the 



2"*"+ wave functions gives the width: 

r(r ^ 77) ^ 0.06 keV . (12) 

6. UNMIXED SCALAR GLUONIA 

6.1. The need for a low mass gb from the sum rules 

Using the mass and coupling of the scalar gluonium G in Table 1 from the USR (Eq.|^), 
into the SSR (Eq.||) sum rules, where 0] ifjgiQ) ~ — 16(/5i/7r)(asG^), one can notice 



27| , p!7| that one needs a low mass resonance ub for a consistency of the two sum rules. 
Using M„g ~ 1 GeV, one gets fg^Jl^]: fag ~ 1 GeV, which is larger than fa ^ .4 GeV. 



6.2. Low-energy theorems (LET) for the couplings to meson pairs 

In order to estimate the couplings of the gluonium to meson pairs, we use some sets of 
low-energy theorems (LET) based on Ward identities for the vertex: 

V{q^ = {p-p'f = 0) ^ {H{pW,\H{p')) ^ 2ml , and V'(0) = 1 , (13) 

and write the vertex in a dispersive form. H can be a Goldstone boson {-K^K^rj^i), a rji- 
f/(l)A-singlet , or a ctb- Then, one obtains the sum rules for the hadronic couplings: 

\ E 9sHHV2fs^2Ml , i J2 9sHHV2fs/Ml^l . (14) 

• Neglecting, to a first approximation the G-contribution, the as and cr^ widths to 
7171, KK, ... (we take ^ 1.37 GeV as an illustration) are |17|: 



T{aB 7T7t) ^ 0.8 GeV , T{a'j^ ^ tttt) ^ 2 GeV , (15) 

which suggests a huge OZI violation and seriously questions the validity of the lattice 



results in the quenched approximation. Similar conclusions have been reached in pq-pO 
For testing the above result, one should evaluate on the lattice, the decay mixing 3-point 
function 1^(0) responsible for such decays using dynamical fermions. 
• The previous LET implies the characteristic gluonium decay (we use Mq ~ 1.5 GeV 
and assume a G- dominance) 



r(G -> r/V) ~ (5 - 10) MeV , li^^M ~ o.22 : gc,, ^ sm9p go,,' ■ (16) 

• Assuming that the G decay into 47r° occurs through aB^B, and using the data for 
/o(1.37) ^ (47r0)5, one obtains 

r(G ^ (Tb(Tb ^ 47r) ^ (60 - 140) MeV . (17) 
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6.3. 77 widths and J /ip ^ '-/S radiative decays 

These widths can be estimated from the quark box or anomaly diagrams 
77 widths of the a, a' and G are much smaller (factor 2 to 5) than r(?7' — 77) ^ 



|. The 

4 keV, 



while B{J/ijj — > 7 0", a' and G) is about 10 times smaller than B(J/ijj 
These are typical values of gluonia widths and production rates ||31||. 



777') ^ 4 10" 



7. UNMIXED SCALAR QUARKONIA 



7.1. The ao(980) 

The ao(980) is the most natural meson candidate associated to the divergence of the 
vector current: dfj,V^{x) = (m^ — md)u{i'^^)d. Previous different sum rule analysis of 
the associated two-point correlator gives p]: M^q ~ 1 GeV and the conservative range 
faQ — (0.5 — 1.6) yieV (/^ = 93 Me^), in agreement with the value 1.8 MeV from a 
hadronic kaon tadpole mass difference approach plus a Oq dominance of the KK form 
factor, and includes the recent sum rule determination [32|. A 3-point function sum rule 
analysis gives the widths p3| , |3^ ,|8[] : 



r(ao Tj-K) ~ 37 MeV 



r(ao ^ 77) ~ (0.3-1.5) keV 



;i8) 



while from SU{d,) symmetry, we expect to have: gaoK+RO — yffl'aojyTr- Analogous sum rule 
analysis in the four-quark scheme Q,^] gives similar values of the masses and hadronic 
couplings but implies a too small value of the 77 width due to the standard QCD loop- 
diagram factor suppressions. The {uu~dd) quark assignement for the ao(980) is supported 
by present data and alternative approaches 



7.2. The isoscalar partner S2 = uu + dd of the ao(980) 

Analogous analysis of the corresponding 2-point correlator gives Ms^ ~ as expected 



from a good SU{2) symmetry, while its widths are estimated to be |T7| , |33 



r(52 ^ TT + TT-) ~ 



120 MeV 



25 

r(^2 ^ 77) ^ 7rr(ao 77) ^ 0.7 keV . (19) 
9 



7.3. The K*{U30) = us and 5*3 = ss states 

The -^"0(1430) is the natural partner of the ao(980), where their mass shift is due to 
SU (3) breakings . An analysis of the S3 over the Kq 2-point functions gives : 



MsJMk; ~ 1.03 ±0.02 



Ss 



1474 MeV 



fs., ^ (43 ± 19) MeV , (20) 



in agreement with the lattice result [^], while the 3-point function leads to 
1(^3 ^ K+K-) ~ (73 ± 27) MeV , TiS^ ^ 77) - 0.4 keV . 



(21) 



In the usual sum rule approach, and in the absence of large violations of the OPE at 
the sum rule stability points, one expects a small mixing between the S2 and S3 mesons 
before the mixing with the gluonium aB- 
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7.4. Radial excitations 

The propreties of the radial excitations cannot be obtained accurately from the sum 
rule approach, as they are part of the QCD continuum which effects are minimized in 
the analysis. However, as a crude approximation and using the sum rule results from the 
well-known channels (p,---), one may expect that the value of can localize approxi- 
mately the position of the first radial excitations. Using this result and some standard 
phenomenological arguments on the estimate of the couplings, one may expect [0: 



Ms/^1.3GeV, r(S'2 ^ vr+vr-) ^ (300± 150) MeV, 1(82 ^ 77) ^ (4 ± 2) keV , 
Ms'^ ^ 1.7 GeV, T{S'^ K+K') ^ (112 ± 50) MeV, T{S'^ 77) ^ (1 ± .5) keV. (22) 
7.5. We conclude that: 

• Unmixed scalar quarkonia ground states are not wide, which excludes the interpretation 
of the low mass broad a for being an ordinary qq state. 

• There can be many states in the region around 1.3 GeV (cr', S3, 6*2), which should mix 
non-trivially in order to give the observed /o(1.37) and /o(l-5) states (see next sections). 

• The /j(1.7) seen to decay mainly into KK if it is confirmed to be a 0+"^ state, can 
be the first radial excitation of the S3 = ss state, but definitely not a pure gluonium. 

8. SCALAR MIXING-OLOGY 

8.1. Mixing below 1 GeV and nature of the (t(1000) and /o(980) 
We consider the two-component mixing scheme of the bare states (ctb, 5*2): 

1/0 >= — sin6's|crB > + cos6's|S'2 > , |(t >= cos^^s|(Tb > -|- sin6's|S'2 > (23) 

A sum rule analysis of the off-diagonal 2-point correlator p7|J^,p^: 

^Uq') ^^Jd'x e''^^0\rp{a,)G\x) ^ m,qq (0)|0), (24) 

u,d,s 

responsible for the mass-shift of the mixed states gives a small mass mixing angle of about 
15°, which has been confirmed by lattice calculations using different unput for the masses 
8] and from the low-energy theorems based on Ward identities of broken scale invariance 



39| , if one uses there the new input values PJT^I of the quark and gluon condensates. 
In order to have more complete discussions on the gluon content of the different states, 
one should also determine the decay mixing angle. In so doing, we use the predictions 
for as, S2 77 obtained in the previous sections and the data r(/o — 77) ~ 0.3 keV. 
Then, we deduce a maximal decay mixing angle and the widths 

l^^l ~ (40 -45)°, 
r(/o TT+TT-) < 134 MeV , 9foK+K-/9fo.+.- ~ 2 , 

r((T -> TT+TT-) ^ (300 - 700) MeV , T{a 77) ^ (0.2 - 0.5) keV . (25) 

The huge coupling of the /o to KK comes from the large mixing with the a. For this 
reason, the /o can have a large singlet component, as also suggested from independent 
analysis [^,^. Extending the previous J/ip 'j + X analysis into the case of the 0, 
one obtains the new result within this scheme ||40|| : 

Br[(p ^ 7 + /o(980)] ^ 1.3 x 10~^ , (26) 
in good agreement with the Novosibirsk data of (1.93 ± 0.46 ± 0.5) x 10~^. 
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8.2. Mixing above 1 GeV and nature of the /o(1.37), /o(1.5), /o(1.6) and /o(1.7) 
As already mentioned previously, this region is quite complicated due to the prolifera- 
tion of states. Many scenarios have been proposed in the literature for trying to interpret 
this region ||^JT7|j28|j29|j38|] . However, one needs to clarify and to confirm the data for 



selecting these different interpretations. We shall give below some selection rules which 
can already eliminate some of these different schemes: 

• The /o(1.37) decay into aa —>■ {4:7i^)s signals mixings with the a, a' and G. 

• The /o(1.5) of Crystal Barrel and the /o(1.6) of GAMS are different objects: 

— Both states like to go into aa and r]'ri, which signals a gluon component. 

— The /o(1.5) coupled to vrvr and KK signals a qq component which may come from the 
S2, S3. Then, it can result from the qq mixings with the a, a' and G, like the /o(1.37). 

— The /o(1.6) couples weakly to vrvr and KK, while the ratio of its rji] and t]'t] is propor- 
tionnal to 1/ sin^^p. These features indicate that it is an almost pure gluonium state, 
which can be identified with the G in Table 1 obtained in the quenched approximation 
(this approximation is expected to better at higher energies using 1/Nc arguments 

• The /o(1.7) is the radial excitation 5*3 of the Ss^ss) state (see section 7). 

9. LIGHT 1 + HYBRIDS 



The experimental situation has been discussed in |44]. The sum rule analysis of the 



spectrum is based on the 2-point correlator ip{q^)H associated to the hybrid currents: 
Ol{x) =: ^^^.A.T.^.G'r : , 0^(x) =: ^7^.A,7.75^,Gr : (27) 

which shows that the lowest state is the l""*". 

• An analysis of the p(l mass and coupling have been done in the past by different 
groups where (unfortunately) the non-trivial QCD expressions were wrong leading 



to some controversial predictions [§]. The final correct QCD expression is given in 



The analysis has been extended recently taking into account the effect of a non-trivial 
1/q^ term in the OPE. One obtains the preliminary results ||48|| : 



^ 1.6 GeV , f~p ^ (25 ~ 50) MeV , {Mp, ^ - Mp ^ 200 MeV , (28) 

where the p' is the radial excitation. One can consider this result as an improvment 
of the available sum rule results ranging from 1.4 to 2.1 GeV, M, though, we cannot 



absolutely exclude the presence of the 1.4 GeV candidate The p'-p mass-splitting is 
much smaller than Mpi — Mp ~ 700 MeV, and can signal a rich population of 1 ^ states 
above 1.6 GeV. The hadronic widths have been computed in |H^,E9|, with the values: 



r(p ^ pvr) ^ 274 MeV , T{p ^ 77r) ^ 3 MeV , 

r(p ^ ri'ir) ^ 3 MeV , T{p nn, KK, r^gr/g) ^ 0{ml) . (29) 



These new |^ and old values of (Mp, fp), indicate that the constraint P5| : 



V'h(O) ^ —asim' , (30) 

is quite inaccurate as it underestimates the value of ipni^) — 2/|M| by a factor 10. An 
independent measurement of this quantity, e.g. on the lattice, can be useful. 
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• One can measure the SU{3) breakings and the mass of the 0(ss) from the difference of 
the ratio of moments, which gives: 



20_2 leOvr^ 

— m„ 

3 ' 9 

a value which is in the range of the lattice results of (1.7 ~ 2.2) GeV 



Ml-M?c^—ml —ms{ss)T^0.3GeV^ =^ ^ 1.7 GeV , (31) 



10. HEAVY 1 + HYBRIDS 



The sum rule predictions for the masses are [p] , p2| : 
M^i,cg) ^ 4.1 GeV , M^(5bg) ^ 10.6 GeV , M^,^^^^^ ^ 6.3 GeV , (32) 
where the two former agree (within the 10% systematics of the two methods) with the 



available lattice values ||53| of about 4.4 GeV [resp. 11 GeV], and are above the i/j-tt [resp. 
T-tt] thresholds. A check of these results is in progress Using the l/M^-expansion 
sum rule, the B*{hug) decays are found to be |]52| : 



T{B* Bin) ^ 250 MeV , r(B* ^ Xvr) ^ 50 MeV . (33) 

X^Bi 

11. CONCLUSIONS 

There are some progresses in the long run study and search for the exotics. Before 
some definite conclusions, one still needs improvments of the present data, and some 
lattice unquenched estimates of the mixing angles and widths which should complement 
the QCD spectral sum rule (QSSR) and low-energy theorem (LET) results. However, one 
can already notice that the simple gg-scheme combined with (large) mixings with gluonia 
(the cr,...) can explain the complex 0"'"+ data. It is a pleasure to thank the organizers of 
Hadron 99 for their invitation to present this review in this exotic country ! 
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